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Mechanisms of spin decoherence are of significance in developing solid state qubits. We performed 
a systematic study on spin decoherence in the compound Ke[Vi 5 As 6042 (D 20 )] ■ 8 D 2 O, using high- 
held Electron Spin Resonance (ESR). By analyzing the anisotropy of resonance linewidths as a 
function of orientation, temperature and held, we study fluctuations in the spin Hamiltonian, and 
Hnd the spin-orbit term as major decoherence source. The demonstrated mechanism can alter the 
lifetime of any spin qubit and we discuss how to mitigate it by sample design and held orientation. 

PACS numbers: 75.50.Xx,71.70.Ej,76.30.-v,03.67.-a 


In solid-state systems, interactions between electronic 
spins and their environment are the limiting factor of 
spin phase lifetime, or decoherence time. Important ad¬ 
vances have been recently realized in demonstrating long- 
lived spin coherence via spin dilution [lj-Q and isolating 
a spin in non-magnetic cages [jj, for instance. Phonons 
are an important source of energy relaxation and thus 
spin-lattice interactions need to be reduced, usually by 
cooling to sufficiently low temperatures. The presence 
of a lattice can be felt by spins through orbital sym¬ 
metries and spin-orbit coupling. An isolated free elec¬ 
tron has a spin angular momentum associated with a 
g-factor g e = 2.00232. In many molecular compounds 
containing 3d elements, the orbital angular momentum 
is quenched and one expects a g-factor only slightly dif¬ 
ferent from g e due to the admixture of excited orbital 
states [U into the ground state caused by spin-orbit cou¬ 
pling. In this Letter, we present observational evidence 
that fluctuations in the spin-orbit interaction can be a 
significant source of spin decoherence. This is accom¬ 
plished by analyzing shape and orientation anisotropy of 
ESR linewidths at 120, 241, and 336 GHz, of the molec¬ 
ular compound K 6 [V l yAsg II 042 (D 2 0 )] • 8D 2 0, V 15 in 
short. This system has been shown to have high co¬ 
herence at low temperatures Igj, has demonstrated co¬ 
herent spin oscillations [H| and interesting out of equilib¬ 
rium spin dynamics due to phonon bottlenecking BEL 
However the details of the spin decoherence are still not 
fully understood. This study points to a new decoher¬ 
ence mechanism, fluctuation in the spin-orbit term, and 
how to optimize the decoherence in spin qubits. 


The V 15 cluster anions form a lattice with trigonal 
symmetry (a = 14.029 A, a = 79.26°, V = 2632 A 3 ) con¬ 
taining two clusters per unit cell. Individual molecules 
have fifteen V lv s = 1/2 ions arranged into three lay¬ 
ers, two non-planar hexagons sandwiching a triangle (see 
Fig. Q3a))- Exchange couplings between the spins in the 
triangle and hexagons exceed 100 K 0, 13] and at low 
temperatures this spin system can be modeled as a tri¬ 


angle of spins 1/2. The Hamiltonian for this three spin 
model is given by EL 

fist = fiz + fij + Hdm (1) 


where TLz represents the Zeeman splitting, Hj is 
the symmetric exchange term, and TLdm is the anti¬ 
symmetric Dzyaloshinsky-Moriya (DM) term (see 14] 
for a detailed formulation), fist eigenvalues are shown 
in Fig. [l]]b) and are used to calculate resonant field posi¬ 
tions B res of the ESR spectra through the method of first 
moments [Tit . Additionally, dipolar interactions between 
molecular spins in the crystal are described by: 


01 _ ..2 m\„ to\o o ~ 3 cos 2 (p pq ) 

fidz — Tb /__j 9p\®)9q\y)SpSq 


( 2 ) 


with g 0 the vacuum permeability, gs the Bohr magne¬ 
ton, 0 the angle between the applied field B 0 and the 2 : 
axis [z is T to triangle plane and is also the symmetry 
c axis of the molecule), d pq the distance between sites 

p and q, g p , q {9) = (g 2 a sin 2 0 + g 2 cos 2 0) , g c , a the g- 

tensor components parallel and perpendicular to the z 
axis, S PiQ the molecular spin, <fi pq the angle between S p 
and d pq . Components of this dipolar Hamiltonian are 
most responsible for fluctuations in the system, as de¬ 
tailed below. 

The linewidth of ESR signals can be significantly 
affected by exchange interactions. In V 15 the intra¬ 
molecular couplings are large and the exchange narrow¬ 
ing effect El collapses the ( 2 / + l ) 15 resonances (/ = 
7/2 for 51 V) into one and it also acts to average out fluc¬ 
tuations related to fist■ This leaves fluctuations in Rdz 
as being the major contributer to spin decoherence. 

There are three possible sources of fluctuation 
in Eq.@, the first being the geometrical factor 
(l — 3 cos 2 (f > pq ) dffq = R pq (t ) since both d pq and (j) pq can 
fluctuate (here, t represents time). This case is described 
by Bloembergen et al. EB (Nuclear Magnetic Resonance 
case) and Kubo and Tomita 0 (ESR case). If R pq (t ) 
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fluctuates randomly, its correlation function decays ex¬ 
ponentially (R(t)R( 0)) = R 2 + r 2 exp(—f/rdip) with a 
Fourier spectrum: 


Jr{v) = r 2 


2 T, 


dip 


1 + 4tt 2 z/ 2 7 


dip 


( 3 ) 


where R is an average value of the geometric term 
Rpqi r is an average size of R(t )'s fluctuations and 
the correlation time T d i P is a characteristic of the random 
motion. This procedure is described generally by Ather¬ 
ton [17j and can be applied to any stationary random 
function that is independent of the time origin. The in¬ 
verse square of the decoherence time T 2 is proportional 
to / J R (v)dv EMI. Therefore, the decoherence rate 
depends directly on r. 

Another fluctuation source comes from thermal exci¬ 
tations to different S z states of S p . q and their effect on 
the second moment of a resonance line has been theoret¬ 
ically analyzed in detail 18 - 23 - Following Kambe and 
Usui flflj . the fluctuations Fourier spectrum is directly 
proportional to a temperature dependent factor: 


KU(T) =< S 2 Z > T ~<S Z > 2 t = S 2 —B s (y) (4) 

ay 


where B s (y ) is the Brillouin function, y = T Z S/T , T z = 
hFmw/kB ( F mw is the microwave excitation frequency), 
and S = 3/2 is the total spin state of the three-spin 
triangle. KU(T) is thus similar to r 2 in Eq. ©) and 
originates from a S 2 (t) correlation function similar to 
< R(t)R( 0) >. This formulation breaks down for tem¬ 
peratures near or below the ordering temperature of the 
spin system. For V 15 this temperature is ~0.01 K EH 
which is well below the temperatures used in this study. 

Fluctuations of g{9) in Eq. (j2j can also reduce the de¬ 
coherence time. Deviation of the ^-factor away from g e 
is due to the spin-orbit interaction (&]: 


g = g e I - 2AA 


( 5 ) 


where g is the g-tensor (diagonal [ g a ,9a,9c ] for V 15 ), I 
is the unit matrix, A is the spin-orbit coupling constant 
and A is a tensor defined in terms of the matrix elements 
of the orbital angular momentum L. In general terms, A 
is the coupling between the ground and excited orbitals 
divided by their energy separation. Here, we focus on 
relative fluctuations in the coupling term with an average 
size £ (assumed isotropic) which induce fluctuations in 
the ^-factor and Rdz■ An average fluctuation in the g- 
factor can thus be written as: 


$9(0) = £ (9(0) ~ 9e) • (6) 

Assuming g{9) is a stationary function with small ran¬ 
dom fluctuations, the temperature dependent correlation 
function of a fluctuating l~Ldz(t) is: 

G dz (t) = (H dz (t)H dz ( 0)) 

= /3(g(t) 9 m 2 KU(T)(R(t)Rm (7) 




FIG. 1. (Color online) (a) Ball-and-stick representation of 
V 15 , with the s = 1/2 V ions in blue. The x axis is along one 
side of the triangle while the z axis is perpendicular to the 
triangle plane and represents the c axis of the crystal unit cell, 
(b) Level diagram of the three spin model with positions of the 
three experimental frequencies shown. Dashed lines show the 
S = 1/2 doublets with the red dashed arrows indicating those 
transitions. Lines show the S = 3/2 quartet with blue arrows 
indicating the transitions; the resonance fields are averaged 
in the first moment calculation of B res at a given frequency. 


where (g(t)g( 0)) = g{9) 2 + (6g(9)) 2 exp(-f/r 9 ), r g is 


the correlation time of ^-factor fluctuations, and /3 = 

/ 87 r 

sity: 


\ 2 . . 

1 ) . As in Eq. ©, this leads to a spectral inten- 
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= T n 


T dip and T gd 


. In ab- 


w / dip ■ 

® reduces to the first 


where V - 's 

sence of g-factor fluctuations, Eq. 

term in the sum as in E3 [see the geometrical term in 
Eq. (|3l) ]. Such a temperature dependence of the linewidth 
has been experimentally demonstrated in Fes 22-24j|. 
nitrogen-vacancy color centers in diamond [ 25 I and other 


studies seem to confirm it as well 2(1 271. If the g -value 
does fluctuate then all three terms in Eq. ® represent 
sources of decoherence. Because T 2 is inversely propor¬ 
tional to / Jdz(v)dv , an important consequence is that 
one can combine different decoherence sources by sum¬ 
ming their effect on J dz (v) as follows: 


n 


E 


rp2 ' 

1 2i 


( 9 ) 


similar to the well-known fact that the sum of uncor¬ 
related variances is equal to the total variance. Addi- 
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tionally, the weight of each term in the sum depends on, 
or can be tuned with, the field angle 0. Here we show 
that for V 15 , fluctuations in the g- value and thus, in the 
spin-orbit interaction, are affecting the coherence time. 

Continuous-wave ESR measurements at 120, 241, and 
336 GHz are performed using the quasioptical super¬ 
heterodyne spectrometer at the National High Magnetic 
Field Laboratory 0,0], with a sweepable 12.5 T super¬ 
conducting magnet (homogeneity of 10~ 5 over 1 cm 3 ). 
Sample temperature can be varied from room tempera¬ 
ture down to 2.5 K. A single crystal of volume < 0.1 mm 3 
was positioned on a rotating stage allowing for continu¬ 
ous change of the angle 9 between B 0 and the c axis 
of the molecule following the procedure described in our 
previous work 0 - The homogeneity of the magnet com¬ 
pared to the size of the crystal allows us to ignore Bq as 
a source of broadening. The applied fields are above 4 T, 
past the crossing of the 5=1/2 doublet and 5 = 3/2 
quartet, such that the ground state of the system is in 
the 5 = 3/2 quartet (see Fig. [U(b)). 

ESR spectra at temperatures T = 4—60 K for B 0 || and 
T to c-axis (9 = 0°,90° respectively) show a Lorentzian 
(homogenous) lineshape. Representative spectra with 
Lorentzian and Gaussian fits are shown in Fig. UK a) 
for comparison. The temperature dependence of the 
linewidth is shown in Fig. H 0 ) for three microwave fre¬ 
quencies Frn.m ■ Compared to measurements made at 
lower fields [30|, where the ground state is in the 5 = 1/2 
doublet, the linewidths are ~ 10 times narrower. Plot¬ 
ted is the full width at half maximum (FWHM) of the 
Lorentzian fits vs T/F mw to underline that the tempera¬ 
ture dependent mechanism of decoherence in the system 
qualitatively follows the temperature behavior predicted 
by KU(T) [Eq. (pT|)]. Additionally, there is no hyper- 
fine structure visible in the spectra (exchange narrowing) 
since the 3d electrons of V interact with the nuclei of 
several other V ions due to the large exchange couplings 
(~ 10 2 K) within the molecule. Due to these properties 
of the measured line width and shape, we can estimate 
Ti to be the inverse of the FWHM. 

There are two distinct curves in Fig. He), dependent 
on the orientation of Bq. To probe this orientation de¬ 
pendence, the linewidth is measured as a function of the 
field angle 9 [see Fig. [21(b)]. The narrowest linewidth 
occurs when 0 = 0° ( Bq || c axis) while the largest oc¬ 
curs at 0 = 90° (Bq in the triangle plane). This im¬ 
plies more decoherence the further g(9) gets from g e since 
g( 0°) = g c « 1.98 and 3(90°) = g a ~ 1.95. This means 
that the FWHM and the first term in Eq.[5]have angular 
dependences that are 180° out of phase. This property 
of the V 15 compound makes it particularly suitable to 
study the effect of spin-orbit fluctuations. Furthermore, 
the fact that the width is largest when g(9) is minimum 
and vice-versa rules out exchange narrowing being the 
cause of this anisotropy since it would require a linewidth 
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FIG. 2. (Color online) (a) Typical measurements of the 
derivative of the absorption x" a t 336 GHz and 241 GHz 
with derivative of Gaussian (blue dotted line) and Lorentzian 
fits (red dashed line), (b) FWHM of Lorentzian fits as a func¬ 
tion of field angle 9 measured at three frequencies: 336 GHz 
(blue triangles), 241 GHz (red circles) and 120 GHz (black 
squares) The dashed lines are calculated widths as a function 
of 9 using Eq. (fill) , (c) FWHM of Lorentzian fits vs temper¬ 
ature/frequency for the 3 studied frequencies. Dashed lines 
are calculated FWHM(T) for 0 = 0° and 90°. 


ex (1 + cos 2 0) [3], as in the case of CsCuCl 3 311. Con¬ 
sequently, this out-of-phase behavior provides strong ev¬ 
idence that 6g(9) ^ 0 and terms in Eq. ([8]) containing it 
must be considered. 

The measured held linewidths can be converted into 
broadening of the molecule eigenvalues through exact di- 
agonalization of the three-spin Hamiltonian, Eq. ([!]). We 
calculate the minimum and maximum excitation frequen¬ 
cies Fmin,max leading to resonance fields Bm ax , min using 
the first moment method 0- Their difference is: 


A E = h(F max - F min ). (10) 

As mentioned above, since I/T 2 oc J Jd z {y) 00 the 
linewidth square is proportional to the prefactor KU(T) 
and either R or r depending on which term in the sum 
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of integrals from Eq. m dominates. This leads to the 
following fit equation [see also Eq. ©]: 


A E 2 


KU(T) A 2 (#)-3 e f + «V(»)(sW 



( 11 ) 

where A and a are fit parameters describing the out- 
of-phase behavior discussed above. Taking AE(0° ) and 
AE(90°) we solve for A and a and calculate the curve 
A E(0) using these values. This curve is then converted 
back to units of Tesla by solving the static Hamilto¬ 
nian. We calculate the resonance fields for two frequen¬ 
cies FE n max = F mw ± AE(6)/2h and take their differ¬ 
ence as the calculated value of the FWHM of an experi¬ 
mental signal. Shown in Fig. [2jb) are the field angle de¬ 
pendencies of the FWHM for the three frequencies used 
in this study as well as calculated widths using Eq. m 

Parameters A and a are then averaged, giving values 
A = 320±34 GHz and a = 10.8±1.6 GHz, and are further 
used to test the validity of model (fill) . The values of A E 
(converted to Tesla) vs T/F mw are fitted with Eq. (fill) 
for two angles 0 = 0° and 90° corresponding to g-factors 
g c = 1.981 and g a = 1.953 respectively. The fit using A 
and a is in very good agreement with the experimental 
data, as shown with dotted line in Fig0c). On the low 
end of the ratio T / F mw one observes a saturation of the 
linewidth, attributed to other decoherence sources. 

Our study provides insight on how to analyze com¬ 
bined sources of decoherence and, in particular, how to 
mitigate the effects of spin-orbit fluctuations. It is evi¬ 
dent from Eqs. © and 0 that the g-tensor should be as 
close as possible to g e . In molecular compounds this can 
be achieved by engineering the ligands type since local 
symmetry affects the diagonal values of the g-tensor of 
a magnetic ion. Aside from material design by chemical 
methods, Jdziy) can be minimized by applying the mag¬ 
netic field at a specific angle 9. For V 15 , this would be 
0 = 0 for which the decoherence times reaches several 
nanoseconds. This time can reach ~ 400 ns by reducing 
r in Jdz{v) via dilution in liquid state, thus allowing the 
observation of Rabi oscillations and spin-echoes [§j. 

In conclusion, we find that the angular dependence of 
the ESR linewidth of V 15 cannot be explained through 
solely a spin flip-flop model [njj], nor through exchange 
narrowing theory |l5| . This dependence is out-of-phase 
with the g-factor orientation anisotropy and we introduce 
a quantitative model based on random fluctuations of 
g. This type of fluctuation can be detrimental to any 
spin-based qubit implementation in solid-state systems 
and we discuss mitigation methods based on chemical 
engineering and experimental parameters. 
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